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IMMERSED FINITE ELEMENT METHOD FOR EIGENVALUE 
PROBLEMS IN ELASTICITY 

SEUNGWOO LEE*, DO Y. KWAK*, AND IMBO SIM t 


Abstract. We consider the approximation of eigenvalue problems for elasticity equations with 
interface. This kind of problems can be efficiently discretized by using immersed finite element 
method (ITEM) based on Crouzeix-Raviart Pl-nonconforming element. The stability and the op¬ 
timal convergence of IFEM for solving eigenvalue problems with interface are proved by adapting 
spectral analysis methods for the classical eigenvalue problem. Numerical experiments demonstrate 
our theoretical results. 
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1. Introduction . In this paper, we consider the approximation of eigenvalue 
problems with interface in elasticity. Eigenvalue analysis is essential basis for many 
types of engineering analysis. As eigenvalues are closely related with the frequency 
and shape of structures, computing the eigensolutions is important to interpret the 
dynamic interaction between the structures. If the frequency of structures is close to 
the system’s natural frequency, mechanical resonance occurs. It may lead to catas¬ 
trophic failure or damage in constructed structures such as bridges, buildings, and 
towers [1]. 

There have been mathematical studies of finite element methods for eigenvalue 
problems. In [2] various computed examples for Laplacian eigenproblems in planar 
regions are studied and there are references to physical problems where the results 
are relevant. For nonconforming approximation of elliptic eigenvalue problems, it is 
shown that the eigenvalues computed by finite element methods give lower bounds 
of the exact eigenvalues whose eigenfunctions are singular in non-convex polygon [3]. 
The guaranteed lower and upper bounds of eigenvalues based on the nonconforming 
finite element approximation are given in [4]. Moreover, let us focus on eigenvalue 
problems in elasticity. A posteriori error estimator for linearized elasticity eigenvalue 
problems is studied in [5]. It is shown that upper and lower estimates for the error of 
eigenpairs are established in terms of a residual estimate and lower-order terms. In 
[6], a method for three-dimensional linear elasticity or shell problems is presented to 
derive computable estimates of the approximation error in eigenvalues. The spectral 
problem for the linear elasticity equations on curved non-convex domains, as well as 
with mixed boundary conditions is considered in [7). Meddahi et al. [8] present an 
analysis for the eigenvalue problem of linear elasticity by means of a mixed variational 
formulation. This method weakly imposes the symmetry of the stress tensor and is 
free from the locking phenomenon. 

When elastic body is occupied by heterogeneous materials, it is known that gov¬ 
erning equations contain the discontinuous material parameters along the interface of 
materials. To simulate such problems by finite element methods, a common strategy 
is to use fitted meshes along the interface. However, this strategy may require a very 
fine mesh near the interface. An alternative approach, proposed in m uni nn nans], 
is an immersed finite element method (IFEM) which can use any meshes indepen- 
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dent of interface geometry. The idea of an ITEM is to construct local basis functions 
to satisfy the interface conditions. For source problems with interface in elasticity, 
Kwak et al. m present a nonconforming IFEM based on the broken Crouzeix-Raviart 
(CR) element [15] . They prove optimal error estimates and provide numerical results 
for compressible and nearly incompressible materials. Computation results of IFEM 
based on the rotated Qi-nonconforming element are reported in m and the related 
work in this direction can be found in m In addition, the spectral analysis of IFEM 
for elliptic eigenvalue problems with an interface is given in m 

In this work, we analyze the spectral approximation of elasticity interface prob¬ 
lems using Pi-nonconforming IFEM and derive the optimal convergence of eigenval¬ 
ues. Moreover, we provide a series of numerical results of the eigenproblems with 
various shapes of interface for compressible and incompressible materials. As a model 
problem, we consider an elasticity eigenvalue problem where the domain is separated 
into two subdomains by interface. The elastic modulus of the material in each subdo¬ 
main is discontinuous along the interface and the eigenfunctions must satisfy certain 
interface conditions. We construct local basis functions to satisfy the jump conditions 
across the interface. Also our local basis functions are based on CR element. It is 
known that CR element does not lock on pure displacement problems [19]. For a trac¬ 
tion boundary problem, the discrete scheme with a stabilization term is introduced 
to overcome locking m- Since interface conditions are related to traction conditions, 
IFEM based on CR elements does not suffer the effects of locking by introducing 
the stabilization term. Furthermore, optimal orders of convergence in the and 
I/^-norms for IFEM are proved in M- Exploiting the ideas of m we formulate the 
discrete scheme with a stabilization term. The proofs for the spectral correctness of 
IFEM are based on the analysis of pn 1^ E51 El]. Introducing a solution opera¬ 
tor, we use spectral properties of compact and self-adjoint operators in Banach space 
[251 [261 Ell EH). In our analysis we adapt the approximation properties of IFEM from 
m to establish the spectral analysis of IFEM. Our proofs for such spectral approx¬ 
imation are very similar to the proofs of [T8| which introduced IFEM to an elliptic 
eigenvalue problem with an interface. 

The outline of this paper is as follows. In the next section, we give a description 
of elasticity eigenvalue problems with interface. In Section 3, we introduce a local 
basis function satisfying interface conditions and formulate an immersed finite ele¬ 
ment method with a stabilization term. Section 4 is devoted to the analysis of the 
spectral approximation which is proved to be spurious-free. In Section 5, we carry out 
numerical experiments for our model problem. The results demonstrate spurious-free 
and locking-free character of IFEM. 

2. Model problem. Let be a connected and convex polygonal domain in M? 
which is divided into two subdomains and by a interface T = fl dQ~ 
(see Figure Erl) . We assume that the subdomains (2+ and are occupied by two 
different elastic materials. Let A and /i denote the Lame coefficients given by 

Eu E 

(1 + z/)(l - 2z/) ’ ^ 2(1 + u) ’ 

where E is the Young’s modulus and u is the Poisson ratio. We note that the coef¬ 
ficients A and fi are 0 < /ii < /i < /i 2 and 0 < A < oo. The constitutive equation 
is related to the displacement field u := {ui) G and the Cauchy stress tensor 
a := {(Jij) G is given by 

t7(u) = 2/ie(u) -j- Atr(e(u))/, 
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where I is the identity matrix of the linearized strain tensor e := (cij) G 

is 


e(u) = i(Vu + Vu^), 
and the usual trace operator tr(€) is 


2 

tr(€) := y] eu. 


For the sake of simplicity, we assume that the density p is a positive piecewise constant 
in subdomains 0+ and From now on, we consider the Lame coefficients A and 
fi as A := A/p and p := p/p. Let us consider the eigenvalue problem for the linear 
elasticity equation with interface, i.e. 


— divt7(u) = in (5 = +,—), 

( 2 . 1 ) 

II 

0 

( 2 . 2 ) 

[<7(u) • n]p = 0, 

(2.3) 

u = 0 on 9^2, 



where and u are the corresponding eigenvalue and eigenfunction, and the symbol 
[•] denotes the jump across the interface F. 



Fig. 2.1. A domain Cl with interface 


We formulate the model problem (EU into the displacement formulation |29| . 
Multiplying v G and applying Green’s identity to model problem (|2.ip in 

each domain we obtain 


/ 


2p €(u) : e(y)dx ■ 


/ 


A div u div w dx — 


/ 

JdCL^ 


cr(u)n • wds = ( 


/ “ 


• vdx, 


2 

e(u) : e(v) = eij(u)ey(v). 
^1 = 1 


where 
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Summing over s = +, — and applying the interface condition (|2.3D . we have the 
following weak formulation 


a(u,v) = w2(u,v), (2.4) 

where 

a(u, v)= / 2/i €(u) : €(v)dx + / Adivudivvdx 

Jq Jq 

and 

Cc;^(u, v) = / U'vdx. 

Jn 

3. Immersed finite element method. In this section, we introduce an im¬ 
mersed finite element method (IFEM) based on Crouzeix-Raviart elements [15]. Let 
{JCh} be the usual quasi-uniform triangulations of the domain ft by the triangles of 
maximum diameter h. Note that an element K G ICh is not necessarily aligned with 
the interface L. For a smooth interface, provided that h is sufficiently small, we are 
able to assume that the interface intersects the edge of an element at no more than 
two points and joins each edge at most once, except possibly it passes through two 
vertices. We may replace T D K hy the line segment joining two intersection points 
on the edges of each K e JCh- We call an element LC G an interface element if the 
interface V passes through the interior of LC, otherwise iF is a non-interface element. 
Additionally we introduce some symbols: 

• - the collection of all interface elements 

• Eh - the collection of all the edges of K e JCh 

We are going to construct local basis functions on each element K of the triangulation 
JCh- For a non-interface element K G JCh, we choose a standard Pi-nonconforming 
basis whose degrees of freedom are determined by average values on each edge of an 
element K. Let N/j,(P) denote the linear space spanned by the six Lagrange basis 
functions 


4^12) 5 J — I5 2, • • • , 6, 

satisfying 

— ^ij -> 

4i2 ds = , 

for each edge Cj of an element P, j = 1,2,3. The Pi-nonconforming space 'Nh{Cl) is 
given by 





^ 0 = (01,02)k ^ ^h{K) for each K e JCh\ JC^; 
if Pi, P 2 ^ JCh share an edge e, then, for i = 1, 2, 

/ 4>i\dK^ds= / and / 0^ ds = 0 

< Je Je JdKndQ > 
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As 



For an interface element K ^ ICh (see Figure [SH]), we describe how to construct 
the basis functions which satisfy the interface conditions (|2.2p . (|2.3p . The piecewise 
linear basis function 0^, i = 1, 2, • • • ,6, of the form 


^iix,y) = 


satisfies 


_ I ^ii 

i2 


(t>i {x, y) = 

0* (a:, y) = ' 


_ / ^0 


Kx- 


■4y 


+ b^x + c'{y 


_ / “0 


box- 


■CoV 


a. + 6-, X + C-, -u 


{x,y) e T+, 


{x,y) e T-, 


I I f bij‘, j 2,3, 

\^j\ Jcj 

1 r [ ^i2 ds 1 j i-5 3, 

lejlJe, 


mm= 0, 

[0i(^)]=O, 


<7($J • n 


= 0 . 


We can express these conditions as a square system of linear equations in twelve 
unknowns for each basis function 0^. It is shown that this system has a unique solution 
regardless of the location of the interface (see [H]). Let us denote ^h{K) as the space 
of functions on an interface element which is generated by i = 1,2,-•• ,6. 
Using this local finite element space, we define the global immersed finite element 
spaee N/i(0) by 


N;,(U) = < 


0 G Nh{K) if T G ICl, and 0 G Nh{K) K ^ /C^; 
if Ki and K 2 share an edge e, then (p = (0i, ^ 2 ) satisfies, 
f^$i\dKids = J^4>i\aK2ds; and Igj^ng^fpids = 0, (i = 1,2) 


In order to describe analysis of IFEM, we introduce some spaces and their norms. 
For a bounded domain D and non-negative integer m, we let = W^{D) be 

the usual Sobolev space of order m with (semi)-norms denoted by || • \\m,D (| * \m,D) 
and let 


:= { u e : u^nn^ € {H^{D n s = +, - }, 
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equipped with norms 


•= I"lm,i3nn+ + |u|m,Dnn-> 
l|u|||„(^) ;= l|u||^^£,nn+ + l|u|lm,Dnn-- 

In addition, we define the space Hhi^) by Hhi^) •= (i^o(r^))^ + 

The ITEM for the eigenvalue problem (|2TD is to find the eigensolution (cj^, Uh) G 
C X Nh{^) such that 

ahiuh,Vh) = \/vh e NhiO.), (3.1) 

where 

a/i(u,v):= 2^ / 2/i e(u) : e(v)(ix + ^ / Adivudivvdx (3.2) 

KeKnKeKn 

eeSh 

The parameter r in the bilinear form a/i(-, •) is a positive constant which is independent 
of the mesh size h. We define the mesh dependent norm || • \\a h oii fbe space H/i(0) 
by 


KeiCh eeSh 

where 

h\\l,K=[ 2fie{v) :e{v)dx+ [ A |divv|^cia;. 

JK JK 

Remark 3.1. The idea of the discrete scheme is motivated from Hansbo and 
Larson JESS- For a source problem without an interface, they prove an optimal con¬ 
vergence of the scheme. For the problem with an interface, Kwak et al. m show the 
scheme yields an optimal result. 

The coerciveness and boundedness of the bilinear form ah{'r) satisfied |14| . 
Theorem 3.1. There exist positive constants Cb and Cc such that 

|a/i(u,v)| < C 6 ||u||a,ft||v||a,ft, VU, V € 
a/t(v,v) > Cc||v ||2 VvgN/,( 11 ). 


4. Spectral approximation. To analyze the spectral approximation, we intro¬ 
duce the solution operator T : (^d(^))^? which associates the solution 

Tf G (i^o (O))^ of the following source problem with every f G 

a(Tf,v) = (f,v), Vv G (ifo(ll))^ 

The operator T is well-defined because unique solvability for every f G is 

shown in [301 ED- It is clear that the operator T is bounded, self-adjoint and compact. 
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In view of the definition of the solution operator T, if (cj^, u) G C\{0} x is an 

eigenpair of (|2.4p . then u) is an eigenpair for the operator T. In a similar way, 

we can define the corresponding discrete solution operator ^ 

by 

ah{Th^,^rh) = (f,v/,), Vv/, G 

with f G Clearly, T/i is also a bounded and self-adjoint operator. Notice that 

an eigenvalue of fbe operator is given by = 1 where is an eigenvalue 
of the discrete problem dsn). 

Before we show the uniform convergence of T/^ to T, we state some assumptions 
which are suggested to analyze the IFEM for the source problem associated with (EU 

in [Hj. 

• (HI). There exists a constant C > 0 such that 

2 M||Tf|| 5 ,(„) + A||divrf||^,(„) < C||f||o,n. 

• (H2). cry(rf) e 

In fact, the hypothesis (HI) implies the regularity estimate which is known when 
the Lame coefficients are continuous on the domain m- On the other hand, such 
estimate for the interface problems is not available to the best of authors’ knowledge. 
The hypothesis (H2) is required to analyze the consistency error of the scheme (j3.2p . 
From now on, we assume the hypotheses (HI) and (H2). 

The following theorem El states the uniform convergence of to T which plays 
an important role in spectral approximation. 

Theorem 4.1. There exists a constant C > 0 such that 

\\Tf - nf\\o,n + h\\Tf - nf\\a,h < Ch^\\f\\o,n, Vf e 


We are going to state the theoretical results of spectral approximation within the 
framework of imiiam. Most proofs of theorems stated below are analogous to [18] 
which deals with the IFEM for elliptic eigenvalue problems. Let us introduce some 
notations for theoretical results. To state the convergence of operators, we introduce 
an operator norm ||T||^(x,v) foi* ^ bounded linear operator L : X ^ Y hy 

Iklk 

The distance between eigenspaces is evaluated by means of distance functions 

disth{x,Y) = inf \\x - y\\a,h, disth{X,Y) = sup dist/i(a:;, T), 

xex,\\x\\a,h=i 

where X and Y are closed subspaces of We denote by cr(T) and p{T) {a{Th) 

and p{Th)) the spectrum and resolvent set of the solution operator T (resp. T/^,), 
respectively. For any z G p(T), the resolvent operator Rz{T) is defined by Rz{T) = 
(z — T)~^ from to or from to and the discrete 

resolvent operator Rz{Th) is defined by Rz{Th) = (z — Th)~^ from and Hh(f^) 

[32|. 

To show that the resolvent operators Rz{T) and Rz{Th) are well-defined and 
bounded, we introduce the following theorem. 



Theorem 4.2. For z G p{T), z 0 and h small enough, there are eonstants 
C'l, C 2 > 0 depending on only O and |^| such that 

||(^ - T)f|U,^ > Ci\\f\\a^, Vf G H^(0), (4.2) 

||(^ - nna,h > C2ma,h. Vf G (4.3) 

Proof. The proof of the first inequality (|4.2p is essentially identical to that of 
Lemma 4.1 from [TB]. The second inequality (|4.3D follows from the first inequality 
(|4.2D and Theorem 14.11 (see Lemma 1 in [22]). □ 

Let ^ be an eigenvalue of T with algebraic multiplicity n and A be a Jordan curve 
in C containing which lies in p{T) and does not enclose any other points of cr(T). 
We define the spectral projection from into (iLo(fl))^ by 

Owing to Theorem 14.21 we can define the discrete spectral projection from 

(1/^(0))^ into H/i(fl) for h small enough by 

We simply denote the projections E{ffj and Eh{C) by E and Eh, respectively. 

Theorem 4.3. The discrete projection operator Eh converges uniformly to the 
projection operator E, i.e., it holds that 

ll-E’ - -B/i||^((L2(n))2,H^(n)) = 0. 

Proof We remark the residual identity 

Rz{T) - R,{Th) = Rz{Th){T - Th)Rz{T), 


so that 

ll-E’ - -B/»||^((L2(n))2,Hh(n)) < 11 -^ 2 ( 1 ) 1)11 jf(Hh(n),H^(n))||r - Th\\^((L^(n)y,nh{n)) 

■ \\Rz{T)\\^((L^Q)y,(L^(Q.))y- 

By Theorem 14.21 and Fredholm alternative [32|, the resolvent operators Rz{Th) and 
Rz{T) are bounded for h small enough. In addition, the operator Th converges to T 
uniformly by Theorem 14.11 Therefore, we conclude the proof. □ 

Finally, we can say that the discrete problem (j3.ip is a spectrally correct approx¬ 
imation of the problem (j2.4p , provided that the following theorem holds m- For the 
proofs of following result, we refer to those of Theorem 1,2,3 and 6 from m 
Theorem 4.4. 

• (Non-pollution of the spectrum) Let A C C be an open set containing cr(T). 
Then for sufficiently small h, cr{Th) C A. 

• (Non-pollution of the eigenspace) 

lim distft(^ft(H^(D)),^((ifi(n))2)) = 0. 
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• (Completeness of the eigenspaee) 

lim disth{E{{Hpn))^),Eh{lihm) = 0. 
h^O 

• (Completeness of the spectrum) For all z G er{T), 

lim d\sih{z,a{Th)) = 0. 

It remains to show the convergence analysis of eigenvalues. The convergence rate 
of eigenvalues is obtained by the spectral properties of compact operators and the 
uniform convergence of the operator to T in Theorem 14.11 

Theorem 4.5. Letf^ he an eigenvalue ofT with multiplicity n. Then for h small 
enough there exist n eigenvalues ofT^ which converge to ^ as follows 

sup \^ - ^i,h\ < Ch"^, 

l<i<n 


where a positive constant C is independent of and h. 

Proof The existence of is a direct consequence of Theorem 14.41 To estimate 
the convergence rate of we introduce some auxiliary operators. Let and T 
be the restriction of operators and T to respectively. Following the 

arguments in [211 [24], we have that the inverse : Eh(H.h{^)) ^((^^(^))^) 

is bounded for h small enough. Hence we can define := and Sh '•= 

Note that the operator Sh is bounded and Sh^ = f for any f G E{{LS{TL)Y). 
The auxiliary operators T, Th^ Sh and ^h provide a following property, for any f G 

E{{L\n)n 

= ShTt - ^pThEht 
= ShTt-^pEhTht 
= Sh{T-n){. 

In view of definition of operator norm dm and Theorem oi we have 


sup 1^ - < C\\T - Th\\^(^E((L^(Q.)Y),E((L'^(Q,))^)) 

l<i<n 


= C sup 

feE((L2(0))2) 

= C sup 

feE((L2(0))2) 

< C sup 

feE((L2(0))2) 

< C/i^ . 


||(T-r;,)f||o,n 

l|f||o,n 

||5^(r-Tft)f||o,n 

llfllo.n 

||(r-Tft)f||o,n 

l|f||o,n 


□ 

Remark 4.1. Overall, we show that our IFEM is spurious-free and has optimal 
convergence property by Theorem j^.^j and Theorem Although uniform conver¬ 
gence of solution operator which is essential basis for spectral analysis is based on the 
hypothesis (HI) and (H2), a variety of numerical results reported in the next section 
corroborate our theoretical results. 
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5. Numerical results. In this section we present a series of numerical experi¬ 
ments to verify the theoretical analysis for the approximation of the eigenvalue prob¬ 
lem (|2.ip in the previous sections. We recall the definition of the Lame coefficients of 
a material 

Eu _ E 

(l + v)(l-2vy ^“2(1 + 1/)’ 

where E is the Young’s modulus and u is the Poisson ratio. We carry out numerical 
tests for the cases of the compressible elastic materials {u < 0.5) and the nearly 
incompressible elastic materials (z/ « 0.5) with various shapes of interface in Figure 
o For a square domain Q = [—1,1]^, we use uniform triangle meshes with mesh size 
h = 2lN where the refinement parameter N is the number of elements on each edge. 
Since analytical expressions for the eigenvalues are not available for all of the examples, 
we use the numerical results on a sufficiently refined mesh as the reference eigenvalues 
in order to estimate the order of convergence. In all the numerical examples, the IFEM 
is implemented in a C++ code and the eigenvalues are computed with ARPACK m- 


D 


ClD 

Circular interface 

Elliptical interface 






© 

Q- 


© © 


Straight-line interface Multiple interfaces 

Fig. 5.1. Domain and interfaces in Examples 1,2,3,4- 5 


Example 1 {Circular interface). In this example, we consider the eigenvalue 
problem (EH with a circular interface. The interface F is a circle with radius r = 0.6 
dividing ft = [—1,1]^ into subdomains and as follows, 

11+ := {(x, y) :x‘^ := {(x, ^) : + / < r^}. (5.1) 

We set the Lame coefficients as (/i“,/i+) = (0.5, 5), (5, 0.5), A+ = 5/i+ and Poisson 
ratio Z/+ « 0.417. Table [5Tl shows the first six eigenvalues and their rates of conver¬ 
gence. The first columns are the reference eigenvalues computed with very fine mesh 
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size h = 2~^ and the other columns are the eigenvalues obtained with IFEM for vary¬ 
ing h. We observe that the convergence rates of the eigenvalue errors are quadratic. 
An eigenfunction for eigenvalue cjI, together with each x and ^-component of eigen¬ 
function, are depicted in Figure [521 


Table 5.1 

First six eigenvalues eomputed by IFEM with eireular interfaee for eompressible materials. 
The referenee eigenvalues in the first eolumn are eomputed with h = 1/2®. The numbers in 

parentheses show eonvergenee rates. 


Circular interface - = (0.5,5), = 5/i^ 


n 

ref 

h = 1/23 

h = 1/2"^ (ord) 

h = 1/25(ord) 

h = 1/2® (ord) 

h = l/2^(ord) 

18.824 

20.409 

19.197 (2.09) 

18.915 (2.03) 

18.846 (2.02) 

18.829 (2.07) 

23.384 

24.040 

23.545 (2.03) 

23.420 (2.16) 

23.392 (2.19) 

23.385 (2.51) 

23.385 

24.953 

23.832 (1.81) 

23.499 (1.97) 

23.413 (2.02) 

23.391 (2.12) 

40.666 

43.609 

41.349 (2.10) 

40.831 (2.05) 

40.706 (2.03) 

40.675 (2.09) 

40.667 

49.381 

42.751 (2.06) 

41.182 (2.01) 

40.795 (2.01) 

40.698 (2.07) 

45.938 

51.811 

47.431 (1.98) 

46.291 (2.08) 

46.023 (2.05) 

45.957 (2.14) 

Circular interface - 

(m" , ) = (5, 0.5), A± = 5^1* 



n 

ref 

h = 1/23 

h= l/24(ord) 

h = 1/2®(ord) 

h = 1/2®(ord) 

h = l/2^(ord) 

7.151 

7.356 

7.205 (1.93) 

7.1652 (2.00) 

7.155 (2.00) 

7.1524 (2.06) 

10.121 

10.165 

10.135 (1.63) 

10.124 (2.21) 

10.121 (1.99) 

10.121 (2.20) 

10.121 

10.313 

10.177 (1.76) 

10.135 (1.98) 

10.124 (1.98) 

10.121 (2.07) 

24.205 

25.696 

24.585 (1.97) 

24.292 (2.12) 

24.224 (2.15) 

24.208 (2.35) 

24.205 

26.708 

24.863 (1.93) 

24.368 (2.01) 

24.245 (2.03) 

24.214 (2.12) 

32.257 

34.569 

32.937 (1.77) 

32.439 (1.90) 

32.303 (1.97) 

32.268 (2.06) 


Example 2 {Elliptical interface). The second example concerns an elliptical 
interface given by T = {(x, y) : x‘^ jo? + = 1} where a = 0.6 and b = 0.3. We set 

subdomain to be an interior and to be the other part of the domain 11, i.e., 

:= {{x,y) : +y'^/b‘^ > 1}, := {(x,y) : x"^/a^ + y"^/b'^ < 1}. (5.2) 

Let Lame coefficients be (/i“, /i+) = (0.5, 5), (5, 0.5) and = 5/i^. In Figure (521 we 
show the errors of the first four eigenvalues by IFEM and corresponding order of con¬ 
vergence. The reference solution is the numerical results on a refined mesh with mesh 
size h = 2“^. Even though the interface becomes shaper than the circular interface, 
the optimal convergence for eigenvalues is obtained. We display an eigenfunction of 
cjI in Figure (521 which is analogous to Figure (521 

Examples {Straight-line interface). We let an interface be a straight line as 
r = {{x^y) : y = 0.5x — 0.2} and the subdomians 11+ and 11“ be 

0+ := {{x,y) : y > 0.5x — 0.2}, 11“ := {{x,y) : y < 0.5x — 0.2}. (5.3) 

Lame coefficients are the same as previous examples, (//“,//+) = (0.5,5), (5,0.5) and 
A+ = 5/i+. In Figure (521 we show the errors of the first four eigenvalues computed 
with IFEM. This figure also presents that the rates of convergence are quadratic. Note 
that in this example the interface meets the boundary of the domain. Nevertheless, 
the order of convergence of the eigenvalue is optimal. Figure (521 shows the computed 
eigenfunction corresponding to uo‘1. 
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Fig. 5.2. The figure above is an eigenfunetion of when {fi = (0.5,5), = 5/i^ in 

Example 1. The figures below are x-eomponent and y-eomponent of eigenfunetion of 




Fig. 5.3. The log-log plots of h versus the relative error of the first four eigenvalues with 
an elliptieal interfaee for the ease of = (0.5,5) (left) and {ii~,= (5,0.5) (right) in 

Example 2. The broken line represents the optimal eonvergenee rate. 
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X 


Fig. 5.4. An eigenfunction of uj^ for the case of an elliptical interface when pi =5, pi^ = 
0.5, = ^pi^ in Example 2. The figures below are x-component of eigenfunction on the left and 

y-component of eigenfunction on the right. 


Example4 {Multiple interfaces) In this case, we solved the problem (|2.ip with 5 
circular interfaces. Let subdomains and be as follows 

Q,~ = U,Li{(a:;, 2 /) : {x - aip + {y - hf < rj, 

= n\n-, 

where (ai,6i) = (0,0), ri = 0.26 and {ui^bi) = (±0.5, ±0.5), = 0.19 for 2 < 

i < 5 (see Fignr dO|) . The Lame coefficients are chosen as follows : (/i“,A“,z^“) = 
(1, 2, 0.33), (/i+, A+, Z/+) = (30, 36, 0.27). Figure [5771 illustrates the error and the rates 
of convergence for cj?, (1 < i < 4) by IFEM. The results in Figure [5771 are in good 
agreement with our theoretical analysis in the previous section. Figure [5^ depicts an 
eigenfunction of uo^. 

Examples {Incompressible materials) To experiment the case of the incom¬ 
pressible elastic materials, we set (/i“,/i+) = (0.5,5), (5,0.5), A^ = 5000/i^ and 
^ 0.4999. We carry out similar numerical experiments with a straight-line in¬ 
terface to demonstrate the locking-free character of our method. The domain Q and 
interface F are the same as Example 3. In Eignre lOl we report the computed errors 
of first four eigenvalues by lEEM. According to Eigure 15.91 it can be seen that the 
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Fig. 5.5. The log-log plots of h versus the relative error of the first four eigenvalues with a 
straight-line interfaee for the ease of (/j,~ ^= (0.5,5) (left) and ) = (5,0.5) (right) in 

Example 3. The broken line represents the optimal eonvergenee rate. 


method has thoroughly locking-free feature for solving the elasticity interface prob¬ 
lems. 
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Fig. 5.6. Eigenfunction of uj‘1 when ^ = (0.5,5), = 5/2^ in Example 3 (above), 

x-component of eigenfunction (below on the left), and y-component of eigenfunction (below on the 
right ). 
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Fig. 5.7. The log-log plots of h versus the relative error of the first four eigenvalues with 
multiple interfaees for the ease of = (1,2)^ = (30,36) in Example 4-- The broken 

line represents the optimal eonvergenee rate. 
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Fig. 5.8. Eigenfunction of uj‘^ when ,A ) = (1,2)^ (/i+,A+) = (30,36) in Example 4 
(above), x-component of eigenfunction (below on the left), and y-component of eigenfunction (below 
on the right). 




Fig. 5.9. The log-log plots of h versus the relative error of the first four eigenvalues with 
incompressible materials (u 0.4999^ for the case of = (0.5,5) (left) and (/i“,/i+) = 

(5,0.5) (right) in Example 5. The broken line represents the optimal convergence rate. 












































































